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With  the  help  of  Domenicali’s  equation  and  the  heat  flux  equation,  a  finite  difference  method  is  directly 
used  to  determine  the  temperature  and  heat  flux  distribution  profiles  in  one  dimension  of  a  thermo¬ 
electric  generator  simultaneously  including  the  Peltier,  Fourier,  Joule,  and  Thomson  effects.  This  calcu- 
lative  method  is  also  used  to  evaluate  the  performance  of  a  thermoelectric  generator,  based  on 
experimental  data  of  thermoelectric  materials  PbTei_yIy  and  PbTe:  Na/Ag2Te.  The  efficiency  of  the 
thermoelectric  generator  as  a  function  of  the  electric  current  is  calculated  and  the  effects  of  structure 
parameters  on  the  performance  of  the  thermoelectric  generator  are  revealed.  The  method  proposed  here 
offers  a  simple  way  to  properly  choose  n-type  and  p-type  thermoelectric  materials  of  a  thermoelectric 
generator. 
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1.  Introduction 

Developing  renewable  energy  and  improving  the  energy  con¬ 
version  efficiency  are  of  important  topics  for  a  sustainable  future 
with  the  development  of  human  activities  [1—3],  Semiconductor 
thermoelectric  generators  (TEGs)  can  convert  an  enormous  amount 
of  unused  waste  heat  to  electricity  and  is  expected  to  play  an 
important  role  in  the  global  sustainable  energy  issue  [4],  Compared 
with  the  conventional  energy  conversion  devices,  thermoelectric 
systems  have  obvious  advantages,  such  as  no  moving  part,  quite, 
compact  and  reliable  for  long  term  supply  [5—7],  TEGs  have  been 
widely  used  in  many  fields.  For  example,  a  solar-driven  TEG,  which 
consists  of  a  TEG  and  a  solar  collector,  has  the  potential  to  partly 
convert  solar  energy  into  electricity  under  terrestrial  operation 
condition  [8,9].  In  Ref.  [10],  a  proof  of  the  principle  experiment  has 
been  designed  to  demonstrate  that  the  solar-driven  TEG  is  a  new 
approach  of  capturing  solar  energy.  Besides,  a  TEG  can  be  integrated 
into  the  bottom  of  a  top  system,  such  as  photovoltaic  modules  [11  ], 
thermophotovoltaic  cells  [12,13],  engine  cycles  [14,15],  and  fuel 
cells  [16],  to  recover  waste  heat,  which  is  another  promising  way 
for  energy  conversion  efficiency  enhancement. 
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The  energy  conversion  efficiency  of  TEGs  is  mainly  determined 
by  the  figure  of  merit  ZT,  which  is  closely  dependent  on  the  Seebeck 
coefficient  S,  electrical  resistivity  p,  thermal  conductivity  k,  and 
defined  as  ZT  =  S2T/(/q>)  [17],  Nowadays,  the  theoretical  and 
experimental  researches  are  mostly  focusing  on  increasing  the 
figure  of  merit  ZT.  Venkatasubramanian  et  al.  reported  that  a  p-type 
Bi2Te3/Sb2Te3  has  ZT  ~  2.4  at  room  temperature  by  controlling  the 
transport  of  phonons  and  electrons  in  the  superlattices  [18],  By 
combination  of  sub-atomic-per-cent  doping  and  nanostructuring, 
Mehta  et  al.  observed  that  bulk  nanomaterials  exhibit  higher  ZT 
than  their  non-nanostructured  bulk  counterparts  [19],  Sevinqli 
showed  through  calculation  that  zigzag  graphene  nanoribbons  are 
very  promising  materials  for  thermoelectric  applications  [20],  It  has 
been  found  from  the  figures  in  Refs.  [21—25]  that  the  Seebeck  co¬ 
efficient  S,  electrical  resistivity  p,  and  thermal  conductivity  k  of 
thermoelectric  materials  are  functions  of  temperature.  It  shows 
clearly  that  the  effects  of  temperature  on  the  material  properties 
should  be  taken  into  account  when  the  performance  of  the  ther¬ 
moelectric  devices  is  evaluated.  However,  only  a  few  researchers 
considered  these  effects.  Domenicali  first  built  the  partial  differ¬ 
ential  equation  through  the  thermodynamic  theory  of  irreversible 
processes  to  find  the  stationary  temperature  distribution  in  an 
electrically  heated  conductor  for  given  boundary  condition  [26]. 
Mahan  pointed  out  that  Domenicali’s  equation  and  a  heat  flux 
equation  must  be  used  to  calculate  the  efficiency  of  a  TEG  when 
transport  coefficients  vary  with  temperature  [27],  This  method 
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allows  the  efficiency  to  be  evaluated  when  a  thermoelectric  device 
is  imposed  any  configuration  of  doping.  Hogan  and  Shih  derived  an 
iterative  technique  based  on  Mahan’s  equations.  A  comparison 
between  the  results  from  the  iterative  technique  and  the  exact 
approach  showed  to  be  in  agreement  with  each  other  [28],  Yama- 
shita  studied  the  effects  of  linear  and  non-linear  components  in 
temperature  dependence  thermoelectric  properties  on  the  cooling 
performance  [29],  The  temperature  dependence  of  material  prop¬ 
erties  were  also  considered  in  solar  TEG  systems.  Li  et  al.  put  for¬ 
ward  a  discrete  numerical  model  to  design  a  concentration  solar 
TEG  [30],  Kraemer  simulated  the  performance  of  a  flat-panel  solar 
TEG,  based  on  the  best  available  properties  of  various  bulk  ther¬ 
moelectric  materials  reported  in  the  literature.  The  best  possible 
conversion  efficiency  of  solar  TEGs  was  calculated  and  was  highly 
consistent  with  the  experimental  measurement  [1],  Xiao  et  al. 
revealed  that  thermal  design  for  a  solar-driven  TEG  could  take  full 
advantages  of  the  characteristics  of  thermoelectric  materials  and 
effectively  improve  the  performance  [31], 

In  this  article,  a  governing  equation  using  finite  difference 
method  for  temperature  in  one  dimension  of  a  TEG  is  originally 
derived  by  employing  Domenicali’s  equation  and  the  heat  flux 
equation.  This  calculative  method  has  quick  convergence  and 
capability  of  searching  the  optimal  solution  within  defined  space. 
Based  on  the  experimental  characteristics  of  the  temperature- 
dependent  Seebeck  coefficient,  electrical  resistivity,  and  thermal 
conductivity  of  samples  PbTei_yIy  and  PbTe:  Na/Ag2Te,  the  tem¬ 
perature  and  heat  flux  distribution  profiles  are  determined 
numerically.  The  efficiency  of  a  TEG  as  a  function  of  the  current 
density  is  calculated  and  the  TEG  is  optimally  designed.  The  effects 
of  structure  parameters  on  the  performance  of  the  TEG  are 
revealed.  When  material  properties  are  independent  of  tempera¬ 
ture,  the  performance  of  the  TEG  can  be  simultaneously  discussed 
by  using  both  the  analytical  method  and  the  finite  difference 
method. 


2.  Governing  equations  for  a  thermoelectric  generator  and 
finite  difference  method 

Fig.  1  shows  the  schematic  diagram  of  a  TEG.  The  length  of  the 
device  is  L.  The  forward  direction  is  defined  from  the  cold  end  to  the 
hot  end.  Tc  is  the  cold  temperature  of  the  device  at  position  x  —  0.  T/, 
is  the  hot  temperature  of  the  device  at  position  x  =  L.JP  and  Jn  are 
the  current  densities  of  the  p-  and  n-type  legs.  The  right  hand  side 
of  Fig.  1  shows  the  control  volume  i  for  a  1  -D  slab,  represented  by  Ax 
at  position  x.  The  interior  nodal  i  is  assumed  to  have  a  constant 
temperature  T)j.  qxi  and  qx,i- 1  are  the  net  heat  flux  densities  into 
nodals  i  and  i-1  for  one  leg,  where  the  subscripts  A  -  n  and  p 
represent  the  n-  and  ptype  materials. 


For  a  thermoelectric  device,  Domenicali’s  equation  is  a  well 
known  equation  to  describe  energy  flux  balance  over  a  segment  of  a 
discretized  thermoelectric  element  at  position  x  [27,28,32,33]  and 
is  given  by 

It  =  -PiWl  gj=;  (A  =  n,p),  (1) 

where  iy(x)  is  the  thermal  conductivity,  T;(x)  is  the  temperature, 
px(.x)  is  the  electrical  resistivity,  and  S;(x)  is  the  Seebeck  coefficient 
along  the  leg.  These  quantities  are  functions  of  position  x  as  the 
temperature  varies  along  the  leg.  The  heat  flux  density  equation  at 
position  x  is  given  by  Refs.  [27,28] 

<?;,(*)  =  kUx)Sx(x)  -  Kx(x)  (A  =n,p),  (2) 

where  the  first  term  indicates  the  Peltier  effect  and  the  second  term 
denotes  the  heat  conduction. 

When  the  Thomson  effect  appearing  in  the  thermoelectric  de¬ 
vice  is  considered,  Eqs.  (1 )  and  (2)  cannot  be  exactly  solved  through 
the  analytic  method.  To  solve  Eqs.  (1)  and  (2),  the  first-order 
approximation  is  adopted,  which  is  similar  to  Refs.  [27,28],  Gov¬ 
erning  equations  for  temperature  should  be  derived  using  energy 
balance  conservation  over  the  control  volume  of  a  half  grid  interval 
width  [34],  The  finite  difference  method  adopted  below  will  be 
proved  to  be  an  effective  method  of  searching  the  optimal  solution. 
Equations  (1)  and  (2)  can  be  converted  into 

- ^  =  c,i  (x)h(x)  +  C;,2(x)q;,(x);  (A  =  n,p)  (3a) 


=  ca(x)Tx(x)  +  cA4(x)q2(x)  +  c25(x);  (A  =  n,p),  (3b) 

where  Qi(x)  =JA(x)/ka(x),  q2(x)  =  -1/*#),  cA3(x)  =  jjSj(x)/KX(x), 
c?a{x)  =  -JxSx{x)Ikx(x),  and  cA5(x)  =  px(x)jj(x). 

By  employing  one  dimensional  finite  difference  method,  Eqs. 
(3a)  and  (3b)  can  be  written  as 

^  Ax'*1^  =  c^.Af  +  c72,i9Ai;P  =  n,P)  (4a) 

and 

^  l  =  C;3-i+1  ^-'+1  +  cHi+i  9«+i  +  c/.5,i+i ;  P  =  n-P) 

(4b) 


Th 


Fig.  1.  The  schematic  diagram  of  a  thermoelectric  generator  and  the  control  volume  for 
a  1-D  slab. 


at  nodal  i.  By  using  Eq.  (4a),  the  heat  flux  densities  into  nodals  i  and 
i  +  1  for  one  leg  can  be  expressed  as 


= 


I-CXU^T  1  T 

C^,fAx  M  cA2iAx  "-1 


;  (A  =  n,p) 


and 


(5a) 


Qii+1  = 


1  -c;ti,-+iAxr 

Ax  Ai+1 


c12,i+lAx 


Txi-,(\=n,p). 


(5b) 


By  using  equations  above,  the  terms  of  heat  flux  densities  can  be 
eliminated  and  the  governing  equation  for  the  temperature  Ty  at 
nodal  i  is  expressed  as 
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aXlTlj-l  +  a>.2^A,i  +  a>3Tu+l  =  =  n-P);  (6) 


where  b)  =  Q5,i+iAx,  a;i  =  1  /c;2.iAx,  a;2  =  (C;4ii+1/%(f+i)- 
(l/C;i2,i+iAx)  -  (1  -  cXyM/cmKx),  and  aB  =  (1  -  cA1_i+1Ax/ 

Q^i+i^)  -cX3,i+rAx  -  (0  -  Qi,i+i^)c;,4:i+i/c/.2,i+i)-  Qm,i(m  = 

1  ~  5)  are  the  material  properties  defined  in  Eq.  (3).  In  this  paper, 
the  Dirichlet  boundary  condition  is  used.  It  means  that  the  tem¬ 
peratures  on  two  sides  of  the  device  are  fixed  [34],  Eq.  (6)  can  be 
applied  to  each  segment  and  iteratively  solved  until  the  desired 
temperature  profile  is  obtained. 

On  the  basis  of  the  determination  of  the  temperature  and  heat 
flux  characteristics,  the  efficiencies  of  the  n-  and  p-type  legs  are 
given  by  Refs  [28,33] 


h  (  jo  S^x)dIj^rdx+J'  j0  PaM d* 


(A  =  n,p) 


pp  =  l/25(Qcm),  and  kp  =  3.194/T(W/cmK)  for  the  p-type  leg  and 
S„  =  0.2687— 329(pV/K),  p„  =  0.1746/(7-310)(Qcm),  and  k„  =  54/ 
7(W/cmK)  for  the  n-type  leg.  The  results  show  that  there  is  no 
difference  between  the  calculated  curves  of  the  present  method 
and  those  of  Ref.  [35],  It  demonstrates  the  validity  of  the  present 
method.  The  main  advantage  of  the  present  method  is  to  simplify 
calculation,  so  that  the  numerical  solution  of  a  governing  equation, 
such  as  Eq.  (6),  can  converge  in  a  matter  of  seconds.  Thus,  the 
present  method  opens  a  more  efficient  way  to  solve  a  similar 
problem. 

In  this  paper,  only  steady-state  behavior  is  studied.  An  impor¬ 
tant  feature  of  the  model  is  its  ability  to  calculate  the  temperature 
profile  taking  the  real  temperature-dependent  properties  of  ther¬ 
moelectric  materials  into  account.  This  feature  is  essential  in 
evaluating  the  performance  of  a  TEG  exposed  to  a  given  tempera¬ 
ture  difference. 

3.  Properties  of  temperature-dependent  thermoelectric 
materials 


and  the  total  efficiency  of  the  TEG  is  given  by  Ref.  [28] 


In  Fig.  2  (a)  and  (b),  efficiencies  for  the  p-type  and  n-type  legs 
are  calculated.  Thermoelectric  properties  can  be  found  by  examples 
from  Chapter  12  in  Ref.  [35],  where  the  Seebeck  coefficient  S, 
electrical  resistivity  p,  and  thermal  conductivity  k  as  a  function  of 
temperature  are,  respectively,  given  by  Sp  —  0.150T  +  211(pV/K), 


(a) 


(b) 


for  given  examples  of  (a)  p-type  material  and  (b)  n-type  material. 


In  this  article,  the  n-  and  p-type  materials,  PbTei_yIy  [22]  and 
PbTe:  Na/Ag2Te  [23],  are  selected  to  evaluate  the  performance  of  a 
TEG,  because  these  materials  have  similar  temperature-dependent 
characteristics  and  relatively  large  figures  of  merit  ZT.  With  the 
coupling  of  large  nanostructures  and  complex  band  structures, 
PbTe:  Na/Ag2Te  becomes  a  high  performance  thermoelectric  ma¬ 
terial.  The  inclusion  of  Ag2Te  in  PbTe  reduces  the  lattice  thermal 
conductivity  by  providing  scattering  phonons  through  nano¬ 
structuring.  The  doping  of  Na  creates  multiple  bands  to  enhance  the 
Seebeck  coefficient  by  improving  electronic  properties  [23],  For 
PbTei_yIy,  Aaron  et  al.  revealed  that  ZT  of  materials  PbTei_yIy  can 
attain  J.4  at  high  temperatures  with  precisely  doping  control  [22], 
Fig.  3  summarizes  the  experimental  characteristics  of  the 
temperature-dependent  Seebeck  coefficient  S,  electrical  resistivity 
p,  and  thermal  conductivity  k  for  different  n-type  samples  PbTei_yIy 
( y  =  0.0004,  0.0007,  0.0012,  0.0020,  0.0028,  0.0040,  0.0055  and 
0.0100)  [22]  and  the  temperature-dependent  Seebeck  coefficient  S, 
electrical  resistivity  p,  and  figure  of  merit  ZT  for  p-type  materials 
PbTe:  Na/Ag2Te  with  different  contents  of  Na  [23],  As  can  be  seen 
from  Fig.  3  (a),  the  absolute  Seebeck  coefficient  and  electrical  re¬ 
sistivity  increase  with  lower  doping  of  1,  while  the  thermal  con¬ 
ductivity  monotonically  decreases.  These  suggest  that  the  doping  of 
I  should  be  optimally  selected  when  PbTei_yIy  is  used  as  a  ther¬ 
moelectric  element.  From  Fig.  3  (b),  it  can  be  found  that  ZT  of  PbTe: 
Na/Ag2Te  increases  with  the  larger  content  of  Na  over  the  entire 
temperature  range,  because  the  electrical  resistivity  is  significantly 
reduced  by  increasing  the  content  of  Na  while  the  Seebeck  coeffi¬ 
cient  and  thermal  conductivity  are  roughly  unchanged.  The  ther¬ 
mal  conductivity  of  PbTe:  Na/Ag2Te  can  be  calculated  from  k  =  S2T/ 
(pZT).  On  the  basis  of  Fig.  3,  the  efficiency  of  the  TEG  will  be  eval¬ 
uated  and  the  TEG  will  be  optimally  designed  in  the  following 
section. 

4.  Performance  evaluation  and  structure  optimum  design 

PbTe  based  thermoelectric  materials  usually  have  high  perfor¬ 
mance  at  middle  temperature  range  ( 600-850  K)  [36],  Because  the 
regions  of  temperatures  given  in  Fig.  3  are  mostly  between  300  K 
and  700  K,  it  is  adopted  in  the  following  discussions  that  the  cold 
end  temperature  Tc  is  300  K  and  hot  end  temperature  is  700  K.  These 
boundary  conditions  are  used  for  both  n-type  and  p-type  materials 
unless  specifically  mentioned.  Based  on  Eq.  (7)  and  experimental 
data  of  n-type  materials  PbTei_yIy  in  Fig.  3  (a),  the  efficiency  t\„  as  a 
function  of  the  current  density  ]n  is  plotted,  as  shown  in  Fig.  4  (a). 
The  length  of  the  leg  is  1  cm.  Fig.  4  (a)  shows  that  the  efficiency  rjn  is 
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Fig.  3.  Experimental  characteristics  of  the  temperature-dependent  Seebeck  coefficient, 
electrical  resistivity,  and  (a)  thermal  conductivity  for  n-type  materials  PbTei  -yly  [22] 
for  different  samples  (x  =  0.0004,  0.0007,  0.0012,  0.0020,  0.0028,  0.0040,  0.0055 
and  0.0100),  and  (b)  figure  of  merit  for  p-type  materials  PbTe:  Na/Ag2Te  [23]  for 
different  contents  of  Na. 


not  a  monotonic  function  of  the  current  density  Jn.  Vn  attains  a 
maximum  value  when  J„  is  adjusted  to  a  reasonable  value.  Gener¬ 
ally,  the  voltage  of  the  thermoelectric  generator  decreases  whereas 
the  current  density  increases,  suggesting  an  optimal  value  for  po¬ 
wer  conversion  efficiency.  Fig.  4  (a)  also  shows  that  the  largest  ef¬ 
ficiency  Vn, max  is  achieved  for  the  sample  with  y  =  0.0007,  because 
this  sample  has  relatively  larger  ZT  over  the  entire  temperature 
range  along  the  thermoelectric  leg.  Fig.  4  (b)  shows  the  efficiency  r}p 
as  a  function  of  the  current  density Jp  for  p-type  materials  PbTe:  Na/ 
Ag2Te  for  different  contents  of  Na.  Jp  is  negative  for  the  forward 
direction  defined  in  Fig.  1.  According  to  Fig.  4  (b),  the  largest  effi¬ 
ciency  r;pimax  is  obtained  for  the  sample  having  the  content,  0.03,  of 
Na,  because  ZT  of  PbTe:  Na/Ag2Te  increases  with  the  increasing 
content  of  Na.  The  samples  PbTe0.9993I0.0007  and  PbTe:  Na/Ag2Te 
with  the  content,  0.03,  of  Na  will  be  used  to  evaluate  the  perfor¬ 
mance  and  design  the  structure  of  a  TEG. 

According  to  Eqs.  (5)  and  (6)  and  the  numerical  method 
mentioned  in  the  second  section,  the  temperature  and  heat  flux  as 
a  function  of  the  position  for  n-  and  p-type  legs  can  be  obtained. 
Fig.  5  shows  the  temperature  versus  position  curves  for  samples 
PbTe0.9993l0.0007  and  PbTe:  Na/Ag2Te.  Fig.  6  shows  the  heat  flux 
versus  position  curves  for  samples  PbTe0.9993I0.0007  and  PbTe:  Na/ 
Ag2Te.  The  current  densities  of  n-  and  p-type  legs  are  Jp  =  -12 A/ 
cm2  and  J„  =  20A/cm2.  It  can  be  seem  from  Fig.  5  that  the 


temperature  profiles  of  two  materials  are  similar  and  not  a  linear 
function  of  position.  It  is  because  the  non-linear  characteristics  in 
the  temperature  dependences  of  thermoelectric  properties  cause 
non-uniform  distribution  of  heat  generation  and  consumption 
along  the  thermoelectric  leg.  Fig.  6  displays  that  the  absolute  value 
of  the  heat  flux  decreases  as  the  position  is  away  from  the  hot  end 
because  more  thermal  energy  is  converted  to  electric  energy. 

According  to  Eq.  (8),  the  efficiency  77  of  the  TEG  can  be  calcu¬ 
lated,  as  shown  in  Fig.  7,  where  the  cross  sectional  area  of  the  p- 
type  leg  is  assumed  to  be  1  cm2  and  the  cross  sectional  area  of  the 
n-type  leg  is  1/3  cm2.  Fig.  7  shows  that  the  efficiency  77  increases  as 
the  current  I  increases,  but  after  reaching  a  maximum  point,  it 
decreases  as  the  current  /  further  increases.  The  reason  is  similar  to 
Fig.  4(a)  and  can  be  explained  by  the  current  voltage  characteristics 
of  a  TEG.  As  the  current  increases,  more  heat  is  converted  to  elec¬ 
tricity  by  the  Peltier  effect,  but  7;  can  not  infinitely  increase  because 
of  the  reduction  of  voltage.  The  maximal  value  ?7max  of  the  efficiency 
is  larger  than  10%,  as  shown  in  Fig.  7.  The  efficiency  can  be  further 
improved  if  the  structure  of  the  TEG  is  optimally  designed. 

According  to  Eq.  (8)  and  Fig.  7,  the  optimization  analysis  is 
applied  to  discuss  the  cross  sectional  area  effects  on  the  perfor¬ 
mance  of  TEG,  as  shown  in  Fig.  8,  where  the  horizontal  coordinate  is 
the  cross  sectional  area  An  of  the  n-type  leg  while  the  cross 
sectional  area  of  the  p-type  leg  equals  1  cm2.  In  Fig.  8,  the  efficiency 
V0pt,A„  has  been  maximized  with  respect  to  the  current.  It  is 
worthwhile  to  point  out  that  VoptAa  is  not  a  monotonically 
increasing  function  of  An.  VoptAn  increases  with  the  increase  of  An 
when  An  is  small.  However,  Vopt^n  attains  its  maximum  (vopt^Jmax 
when  A„  —  A„,m.  If  A„  further  increases,  v0pt,An  will  decrease. 


Jp(A/cm2) 


Fig.  4.  Efficiency  versus  current  density  curves  for  (a)  n-type  materials  PbTei_yIy  and 
(b)  p-type  materials  PbTe:  Na/Ag2Te. 


S.  Su  et  al  /  Energy  70  (2014)  79-85 


-7.6 

"g  -8.0 

t-8.2 


0.0  0.2  0.4  0.6  0.8  1.0 


t 

-6.6  ^ 


x(cm) 


Fig.  6.  Heat  flux  versus  position  curves  for  samples  PbTe0.9993I0.0007  and  PbTe:  Na / 
Ag2Te. 


5.  Two  calculative  methods  of  a  simplified  model 


By  using  the  boundary  condition  1/(0)  =  Tc  andT/(L)  =  77,,  Eq.  (9) 
can  be  integrated  to  obtain 

= “p^(x  - i/2) + K*(Yr*W =  n,p)-  (10) 

Using  Eqs.  (2)  and  (10),  one  can  derive  the  heat  flow  densities  at 
x  =  0  and  x  =  L  as 

<J/(0)  =  y/TcS/  -  p)jjl/2  +  K^Th^Tc).  (X  =  n,p)  (11) 

and 

91(1)  =  hJhSx  +  pjjL/2  +  K,(Th~  Tc);  (X  =  n,p).  (12) 

It  can  be  derived  from  Eqs.  ( 1 1 )  and  ( 1 2 )  that  the  total  heat  flows 
from  the  hot  end  to  the  TEG  and  from  the  TEG  to  the  cold  end  are, 
respectively,  given  by 


In  order  to  conveniently  determine  the  performance  of  a  TEG,  a 
simplified  method  is  often  used  [35,37,38]  by  assuming  that  the 
Seebeck  coefficient,  electrical  resistivity,  and  thermal  conductivity 
of  the  TEG  are  independent  of  temperature.  In  such  a  case,  the 
Thomson  coefficient  equals  zero,  i.e.,  T(x)dS;(x)/dT  =  0;(A  =  n,p).  Eq. 
(1 )  can  be  simplified  into  [35] 

^=-P^  =  n,p).  (9) 


qh  =  SIaTh  - 1 12aR  +  K(Th  -  Tc)  (13) 

and 

<?c  =  SIaTc  + 1/2 R  +  K(Th  -  Tc),  (14) 

where  S  =  (SP-S„),R  =  ( ppL/Ap  +  p„L/An),  and  K  =  ( KpAp/L  +  K„An/L) 
are  the  Seebeck  coefficient  S,  electrical  resistance  R,  and  thermal 
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Fig.  7.  The  efficiency  ij  versus  current  I  curve  of  the  thermoelectric  generator  con¬ 
sisting  of  PbTe0.9993I0.0007  and  PbTe:  Na/Ag2Te. 
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conductance  K  of  a  thermoelectric  device,  and  la  —  Jxflx  is  the 
electric  current  of  a  TEG.  By  using  Eqs.  (13)  and  (14),  the  power 
output  Pa  and  efficiency  rja  of  the  TEG  can  be  analytically  expressed 
as 


Pa  =  qh-qc  =  SIa(Th-Tc)-I2aR 


(15) 


and 


Pa  =  SIa(Th  -  Tc)  -  IgR 
Qh  SIaTh  -  \l2a R  +  K(Th  -  Tc) ' 


(16) 


Fig.  9  shows  the  efficiency  p  as  a  function  of  the  power  output  P 
obtained  from  Eqs.  (15)  and  (16),  where  the  dots  on  the  curve  are 
obtained  by  the  finite  difference  method.  In  Fig.  9,  the  cross 
sectional  area  of  the  p-type  leg  is  1  cm2,  the  cross  sectional  area  of 
the  n-type  leg  is  0.5  cm2,  and  the  average  values  of  the  experi¬ 
mental  data  of  samples  PbTe0.9993I0.0007  and  PbTe:  Na/Ag2Te  ob¬ 
tained  from  Fig.  3,  which  are  listed  in  Table  1,  are  adopted.  It  is 
seen  from  Fig.  9  that  the  results  obtained  by  two  calculative 
methods  are  identical.  The  calculations  of  error  analysis  by  the 
finite  difference  method  and  the  analytical  method  are  also  shown 
in  Table  2.  The  error  is  quantified  as  a  percent  discrepancy,  which 
is  the  ratio  of  the  efficiency  difference  between  the  finite  differ¬ 
ence  method  and  the  analytical  method  to  the  efficiency  calcu¬ 
lated  by  finite  difference  method.  The  largest  error  is  only  1.68%  at 
current  25A.  In  addition,  Fig.  9  shows  that  the  t)  ~  P  curve  is  a 
closed  loop  line  passing  through  the  zeros  point.  On  the  curve, 
there  are  a  maximum  efficiency  Tjmax,  whose  corresponding  power 
output  is  Pm,  and  a  maximum  power  output  Pmax  whose  corre¬ 
sponding  efficiency  is  rjm.  Thus,  the  ranges  of  the  efficiency  and  the 
power  output  must  be  constrained  by  r/m  <  rj  <  r/max  and 
Pm  <  P  <  Pmax-  Obviously,  rjmax  and  Pmax  determine  the  upper 
bounds  of  the  efficiency  and  power  output,  while  rjm  and  Pm  give 
the  lower  bounds  of  the  optimized  efficiency  and  power  output. 
Thus,  rj  max.  Pmax.  Vm,  and  Pm  are  some  four  important  parameters 
ofTF.Gs. 


6.  Conclusions 

With  the  help  of  Domenicali’s  equation  and  heat  flux  equation, 
a  finite  difference  method  is  directly  used  to  determine  the  tem¬ 
perature  and  heat  flux  distribution  profiles  of  a  thermoelectric 
generator  in  one  dimension.  Efficiencies  of  the  thermoelectric 
generator  under  different  conditions  are  calculated  based  on 
experimental  characteristics  of  the  temperature-dependent  See- 
beck  coefficient,  electrical  resistivity,  and  thermal  conductivity  for 
samples  PbTei_yIy  and  PbTe:  Na/Ag2Te.  The  efficiency  and  struc¬ 
ture  parameter  are  optimally  designed.  It  is  proved  that  the  finite 
difference  method  adopted  here  is  accurate  and  quick  conver¬ 
gence  to  evaluate  the  performance  of  the  TEG  by  comparing  with 
other  methods.  The  method  offers  a  simple  way  to  properly 
choose  thermoelectric  elements  of  a  thermoelectric  generator  on 
the  basis  of  the  experimental  characteristics  of  thermoelectric 
properties. 


The  average  values  of  temperature-dependent  material  parameters  for  samples 
PbTe0.9993io.ooo7  and  PbTe:  Na/Ag2Te. 


Sample  Type  a  (gV/K)  p  (m£2  cm)  k  (W/m-K) 


PbTe0.9993l0.0007  n  -189  1.67  1.53 

PbTe:Na/Ag2Te  p  253  3.72  0.947 


Table  2 

The  calculations  of  error  analysis  by  the  finite  difference  method  and  the  analytical 
method. 


Current  (A)  Finite  method  Analytical  method  Error  (%) 


Efficiency 

1  0.0237 

4  0.0738 

10  0.110 

12  0.109 

17  0.0862 

19  0.0699 

22  0.0390 

25  0.000595 


Efficiency 

0.0237  0 

0.0739  0.136 

0.111  0.909 

0.110  0.917 

0.0870  0.928 

0.0707  1.14 

0.0395  1.28 

0.000605  1.68 


Acknowledgments 

This  work  has  been  supported  by  the  National  Natural  Science 
Foundation  (No.  11175148),  973  Program  (No.  2012CB619301),  and 
the  Fundamental  Research  Fund  for  the  Central  Universities  (No. 
201312G007),  People’s  Republic  of  China. 

Nomenclature 

S  Seebeck  coefficient  (pV  K_1) 

ZT  figure  of  merit 

x  position  along  TEG 

L  length  of  TEG 

Tc  cold  end  temperature 

Th  hot  end  temperature 

Jp  current  density  of  p-type  leg 

Jn  current  density  of  n-type  leg 

Ax  control  volume 

Tn  temperature  at  nodal  i 

q;j  heat  flux  density  into  nodal  i 

Tj(x)  temperature  at  x 

q;,(x)  heat  flux  density  at  x 

Qm,j(m  =  1  ~  5)  material  properties  defined  in  Eq.  (3) 

QAi.a72.Q23  factor  of  the  matrix  defined  in  Eq.  (6) 
bx  factor  of  the  matrix  defined  in  Eq.  (6) 

A„  cross  sectional  area  of  n  type  leg 

R  electrical  resistance 

K  thermal  conductance 

I  electric  current 

Pm  power  at  maximum  efficiency 

Creek  symbols 

p  electrical  resistivity  (mfi  cm) 

k  thermal  conductivity  (Wm_1  K-1) 

rjx  efficiencies  of  n-  or  p-type  legs 

i)  total  efficiency 

r)m  efficiency  at  maximum  power 

Superscripts  and  subscripts 

h  hot  end  of  TEG 

c  cold  end  of  TEG 

i  interior  nodal  along  TEG 

n  n-type 

p  p-type 

X  n  or  p  type  material 

max  maximum  value 

opt  optimal  value 

y  different  samples  of  PbTei_yIy 

a  a  parameter  independent  of  temperature 
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